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A nonlinear model of modulational processes in the subsonic regime involving a lin- 
early unstable wave and two linearly damped waves with different damping rates in 
a plasma is studied numerically. We compute the maximum Lyapunov exponent as a 
function of the damping rates in a two-parameter space, and identify shrimp-shaped 
self-similar structures in the parameter space. By varying the damping rate of the 
low-frequency wave, we construct bifurcation diagrams and focus on a saddle-node bi- 
furcation and an interior crisis associated with a periodic window. We detect chaotic 
saddles and their stable and unstable manifolds, and demonstrate how the connec- 
tion between two chaotic saddles via coupling unstable periodic orbits can result in 
a crisis-induced intermittency. The relevance of this work for the understanding of 
modulational processes observed in plasmas and fluids is discussed. 
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I. INTRODUCTION 



The theory of nonhnear modulational processes can explain the generation of modu- 
lated wave envelopes commonly observed in laboratory plasmaJ^, space and astrophysical 
plasmaJ^and oceans^. For example, modulational Langmuir processes in plasmas refer to 
the mechanism by which initially small fluctuations of electrostatic Langmuir waves produce 
density accumulation towards regions where the electric field is weaker. The high-frequency 
Langmuir waves become trapped in cavitons where the plasma density is lower, resulting in 
the formation of Langmuir envelope solitons. If the phase velocity of the density perturba- 
tion is lower than the ion-acoustic velocity, this process is known as a subsonic modulational 
proces^. Modulational processes have been observed in laboratory experiments of non- 
magnetized plasmas in the presence of an external electric fielcP and in Nd-laser beat-wave 
experiments in a homogeneous plasmeP. In space plasmas, modulational processes can ex- 
plain the detection of modulated wave envelopes observed in the auroral regiorP^. 

Theoretical models of modulational processes in astrophysical plasmas have been derived 
and studied numerically by several authors. For example, Pataraya and Melikidze^ used 
the nonlinear Schrodinger (NLS) equation with nonlinear Landau damping to study the 
modulational instability of nonlinear waves in relativistic plasmas. Weatherall'^ studied a 
nonlinear model of plasma waves in pulsars, and found that the onset of wave turbulence 
occurs via a modulational instability. Robinson et al.'^ proposed a unifying framework 
which incorporates modulational and decay processes driven by a pump Langmuir wave. 
They performed analytical and numerical studies to obtain quantities such as growth rates, 
stability boundaries and instability thresholds for various types of modulational and decay 
processes. Recently, Asenjo et al.'^ derived a NLS equation to describe the self-modulation 
of nonlinear waves in a hot magnetized relativistic electron-positron plasma. 

In the weakly nonlinear regime, modulational processes in plasmas and fluids can be 
described by low-dimensional models of nonlinear wave-wave interactions. A great amount 
of literature has been dedicated to their study, for example, Russell and Ott'^l performed 
a numerical study of a three- wave truncation derived from the NLS. By varying a control 
parameter, they found a variety of dynamical behaviors such as transitions from stable sta- 
tionary points to periodic orbits via Hopf bifurcations, period-doubling bifurcations leading 
to chaos, and transitions from chaotic to periodic solutions via tangent bifurcations. Gosh 
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and Papadopoulos'^ showed that the model of Russell and Ott'^l can also apply to a three- 
wave truncation of the derivative nonlinear Schrodinger equation (DNLS) to study the onset 
of Alfvenic turbulence. DuBois et al.'^ performed Monte Carlo simulations of a truncated 
model to study the transition from the sonic to subsonic regime of the modulational insta- 
bility. Chian and Abalde'^ derived a nonlinear model of the hybrid stimulated modulational 
instability where Langmuir, electromagnetic and ion-acoustic waves are coupled, and its tem- 
poral dynamics was studied by varying several parameters such as the frequency mismatch, 
dissipation rate and wave dispersion. In Miranda et al.^ two different types of intermit- 
tency were found in chaotic time series obtained from a three-wave truncation of the NLS. 
Chaos in a three-wave truncation model of the derivative nonlinear Schrodinger equation 
for nonlinear Alfven waves has been studied by Sanmartin et al.^ and Sanchez- Arriaga et 
a} [T71|i8P9 | Mori and Janssen^^ derived an expression relating the kurtosis of the sea surface 
elevation and the occurrence probability of freak waves generated via nonlinear four-wave 
modulational interactions. Their statistical theory is in good agreement with observational 
data of extreme wave eventd^^'^. 

In this paper we study a low-dimensional model of nonlinear four-wave subsonic mod- 
ulational processes. We start our numerical analysis by exploring a two-dimensional (2D) 
parameter space in terms of the maximum Lyapunov exponent. A number of papers have 
shown the existence of regular, self-similar structures in 2D parameter space corresponding 
to periodic oscillations, embedded within regions corresponding to chaotic dynamics. These 
structures are called "shrimps" due to their characteristic shape consisting of a main cen- 
tral body and four elongated "antennae" or "legs"!^^!!!^ g^^^j g^j^g globally organized along 
straight lines, ellipses, and spirals. A theoretical explanation of the spiral organization of 
shrimps was presented by Stoop et al.^^. They performed numerical simulations and lab- 
oratory experiments of the Nishio-Inaba electronic circuits to demonstrate that the spiral 
emergence of shrimps is due to the presence of a homoclinic saddle-focus point in the pa- 
rameter space. Shrimps have been found in several models including discrete-time mapd^^'^, 
continuous-time models represented by ordinary differential equationd^^'^ and in a number 
of experiments using nonlinear circuit^^^*^. Our analysis is then focused on a periodic win- 
dow, in which we perform a numerical detection of chaotic saddles, which are nonattracting 
chaotic sets responsible for chaotic transients in nonlinear systemd^^^. They also play an 
important role in the intermittent dynamics of low- and high-dimensional systems as demon- 
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strated by recent papersS^HMl^ hence chaotic saddles are of fundamental importance for the 
characterization of the system dynamics. 

This paper is organized as follows. In section |Tl] we review the derivation of the nonlinear 



model used in this work. In section III we describe the main numerical algorithms used for 



the computation of the Lyapunov exponents, the detection of chaotic saddles and their stable 



and unstable manifolds. The numerical analysis is then presented in section |lVj We begin by 
inspecting the parameter space and identifying a periodic window in the bifurcation diagram, 
then focusing on a saddle-node bifurcation and an interior crisis. Finally, the conclusion and 
the relevance of our results for the study of modulational processes in plasmas and fluids is 
presented in section |Vj 

II. NONLINEAR MODEL 

The nonlinear interaction between high-frequency Langmuir waves and low-frequency 
ion-acoustic waves in plasmas is described by the Zakharov equationd^ 



^ + u,l,E, - 7.<V2Ej = -i^l^^Ej, (1) 

where E/j represents the high-frequency electric field fluctuations, Upe is the electron plasma 
frequency, 7e is the ratio of the specific heats for electrons, Vth is the electron thermal velocity, 
hq is the average ion density, rii represents the ion density fluctuations (i.e., the deviations of 
the ion density from no), is the ion-acoustic velocity, eo is the permittivity of vacuum, rrii 
is the ion mass, and the angular brackets represent the fast-scale time average (i.e., the time 
average performed in a period corresponding to the high-frequency fluctuations). Assuming 
that the Langmuir waves propagate along the x axis we can write the electric fleld in the 
amplitude-modulated representation as 

Eh{x,t) = ^E{x,t) exp{-iujpet) + c.c., 

where E represents a slowly-varying complex envelope such that \d^E\ ^ \ijjpedtE\, \dlE\ <^ 
\kdxE\, and i = a/— 1. By deflning the following dimensionless variables 
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where rj = {jeKsTe + ''yiKsTi) / {KBTe), Kb represents Boltzmann's constant, and Tj are 
the electron and ion temperatures respectively, 7^ the ratio of specific heats for ions, me the 
electron mass, and Ad represents the Debye length, we can write Eqs. as 



.dE d^E 

dr dz"^ 

d^n d^n d\E\'^ 



^^ + ^ = ^E, (3) 

1712 

(4) 



Q^2 Q^2 

Next, we make the approximation dn/dr = in Eq. Q, which can be justified as 
follows. Random fiuctuations of the high-frequency electric field will produce undulations 
or "ripples" of Langmuir waves. These ripples imply a small gradient of wave intensity 
generating a force called the ponderomotive force, which moves electrons and ions towards 
regions where the intensity of the electric field is minimum. As a result, the plasma density 
will increase within these wave troughs. The dispersion relation of Langmuir waves 
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indicates that, for a given frequency ul, large ki wavenumbers can only exist in regions 
where Upe is small. Hence, Langmuir waves will be "trapped" in low-density regions. This 
localization of Langmuir waves allows the wave envelope to develop a growing ripple, which 
enhances the effect of the ponderomotive forc^^. In the frame of reference of the high- 
frequency Langmuir wave, the low-frequency density perturbation becomes practically static. 
Hence, we can assume that the time derivative in Eq. ^ vanishes. 

Introducing the static approximation dn/dr = in Eq. Q, integrating twice and setting 
the constants of integration equal to zero yield 
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n = -\E\\ (5) 
By substituting Eq. ([S]) into Eq. ^ we obtain 

f ^S + l-P^-O,. (6) 

Equation ([6]) is called the Nonlinear Schrodinger equation because it resembles the Schrodinger 
equation of quantum mechanics^. Following Russell and Ott^^ we include the term ivE in 
Eq. ([6]), where z> represents the linear growth/damping rate of and obtain 

fdE _\ d'^E _,2. 



We look for solutions of Eq. ([T]) in the form of three travelling waves 



E{z, t) = Eq{t) exp[-i(A;o2 - ujqt)] 
+Ei{t) ey.\;)[—i{kiz — Uit)] 

+E2{T)e-xp[-i{k2Z - U2t)], (8) 

where the resonant condition 2kQ = ki + k2 and the linear dispersion u^j = {a = 0, 1, 2) are 
assumed. Introducing (|8]) in ([T]) we obtain the following set of complex ordinary differential 
equations (ODEs) 



^0 = -Kko)Eo + ^[l^or^o + 2|Ei|2Eo 
+2|EopEo + 2EqEiE2 exp(2i(5r)], 

El = -v{ki)Ei + i[\Ei\''Ei + 2|Eo|'Ei 
+2|E2|^Ei + E2E2exp(-2z(5r)], 

E2 = -0{k2)E2 + i[\E2\^E2 + 2|EopE2 
+2|EipE2 + ^o^^i exp(-2z5r)], 



(9) 
(10) 

(11) 



where 5 = {Si+S2)/2, = Wo-— (c^ = 1, 2) and Q = d/dr. In terms of quantum mechanics, 
Eqs. (|9])-(11) denote the nonlinear temporal evolution of a 4-quanta system wherein a pair 
of pump quanta interact with a pair of Stokes and anti-Stokes daughter quanta. This set of 



equations can be rewritten using the polar representation = a^it) exp[i'?/'o-(t)] (cr = 0, 1, 2), 
where a and ip are real amplitude and phase respectively, to obtain a set of four real ODEs 



clq = I'oao + 2aoa.ia.2 sine/), 



Oi = — z/itti — aQa2Sini 



02 



-U2a2 — aQtti sm (p, 
-26 + al + al- 2al 



+ 



4aia2 — Oq 



— + — 

Oi 02 



COS( 



(12) 
(13) 
(14) 

(15) 



where (pit) = 2ipQ — i'l ~ ip2 — 25/!:, Uq = —uiko), 2 = i'{ki^2)- Wave represents a linearly 
unstable (pump) wave, and waves 1 and 2 represent linearly damped (daughter) waves. 
By assuming equally damping rates for the daughter waves {v2 = vi = v) one can obtain 



oi = 02, reducing Eqs. (12)-(15) to a set of three ODEs 



do = t^o^o + 2aQai sin 0, 
di = — z/fli — a^ai sin0, 
= -25 + 2{al - al) + 2(2a^ - a^) cos. 



(16) 
(17) 
(18) 



which has been extensively studiecP^HMSl paper we avoid this simplification and let 



^ ^2 ^ retaining the four-dimensional set of equations (12)-(15). Numerical integration of 



Eqs. ( 12 )-( 15 ) is performed using the Isoda package which is a variable-step integratoi'^. In 



this paper we define a Poincare surface of section as 02 = 1, and refer to "Poincare points" 
as the intersection between the trajectory and the Poincare surface with d2 > 0. 

III. NUMERICAL ALGORITHMS 



A. Computing Lyapunov exponents 

The dynamical stability of a nonlinear system represented by a set of ODEs can be 
characterized by the Lyapunov exponents, which represent the exponential divergence or 
convergence rates of nearby trajectorie^^^*^^. Let us rewrite Eqs. (12)-(15) in the form 
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du 

~dt 



Fiu) 



where u = (oq, cti, 02, 0). Denote by uq an initial condition of system (19) at t 
6'i(uo,to) its solution, that is 



(19) 

to, and 



Oti^o, to) 



iuo,to),t), 



(20) 



and 6tg{uo,to) = uq. Taking the derivative of Eq. (20) with respect to uq and defining the 
flux Jacobian 6i(uo,to) = -Du„6't(uo, to) we obtain 



et(uo,to) = D^F{9t{uo,to),t)et{uo,to), 



(21) 



and 0tQ(uo,to) = ^, where / represents the identity matrix. Equation (21) is known as the 
variational equatiorP^SSI \ small perturbation 6uo of uq will evolve as 



6u = etiuo,to)5uo. (22) 
Suppose that 5u can be decomposed into a linear combination of j orthogonal vectors 



6u = 6ui + 6u2 + ... + Su 
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1 



(23) 



where D is the dimension of the phase space, and 5\io = 6ui{to) + 5u2(to) + ••• + Suj{to). 
The asymptotic behavior of perturbation 6u is given by the Lyapunov spectrum, which is 
the set of Lyapunov exponents Xj defined bj^ 



Xj = lim 



t->oo t — tn 



In 



rrij (t) 



mj{to) 



j = l,...,D. 



(24) 



where rrij are the eigenvalues of the flux Jacobian G^, obtained by solving Eq. (21). A pos- 



itive Lyapunov exponent indicates that the distance between two trajectories will increase 
exponentially in time, whereas a negative value indicates that the distance will decrease 
exponentially in time. The number of Lyapunov exponents is equal to the dimension of the 
phase space, and each exponent represents the exponential divergence/convergence rate of 
initial conditions along a different direction in the phase space. The maximum Lyapunov 
exponent represents the rate in the direction of maximum divergence, or minimum conver- 



gence. Note that, in the case of continuous-time nonlinear models such as Eqs. (12)-(15), 



at least one exponent is always equal to zero and represents the dynamics in the direction 
tangent to the flow, except when the asymptotic solutions converge to a fixed point. In this 
paper we will refer to Amax as the maximum non-zero Lyapunov exponent. Chaos occurs 
when Amax > 0. If Amax < 0, then the dynamics can be periodic or quasiperiodic. Periodic 
behavior occurs when the Lyapunov spectrum is characterized by only one vanishing Lya- 
punov exponent and Amax < 0, whereas quasiperiodic behavior occurs when there are two 
or more vanishing exponents and Amax < 0. 

The numerical techniques for computing the Lyapunov spectrum of ODEs were first 
developed by Shimada and Nagashima'^ and Benettin et al.'^ . In this paper we compute 



the Lyapunov spectrum of Equations ( 12 )-( 15 ) following the numerical implementation found 
in Appendix A of Wolf et al.*^. 



B. Detecting chaotic saddles 

Since chaotic saddles are non-attracting chaotic sets, the trajectory of any initial condi- 
tion arbitrarily close to the chaotic saddle will eventually leave its vicinity in forward-time 
dynamics, except for initial conditions located exactly in the chaotic saddle or its stable 
manifold. The stable manifold is defined as the set of points which converge to the chaotic 
saddle in forward-time dynamics, whereas the unstable manifold is the set of points which 
converge to the chaotic saddle in reverse-time dynamic^^^'^. We implement the following 
two methods to numerically find a chaotic saddle. 



1. The PIM-triple method 

This method, proposed by Nusse and Yorke'^, defines a "Proper Interior Maximum 
triple" , or PIM-triple, as the triple (a, c, b) which lies on the line segment (a, b) intersecting 
(or "straddling" ) the stable manifold of a chaotic saddle in phase space in such a way that 
point c is closer to the stable manifold than a and b, and the escape time profile of the triple 
(a, c, b) will have a maximum at c. This procedure can be described as follows 

1. Define a region of the phase space with no attractors, called the restraining region. 
Define the escape time as the time it takes for a trajectory to leave the restraining 
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region. Select two points a and b on the Poincare surface which straddle the stable 
manifold of the chaotic saddle in phase space, and set the error e to a small value. 



2. Divide the segment (a, b) into M subsegments, and compute the escape time of each 
subsegment. Choose one subsegment (ai,ci,6i) which qualifies as a PIM-triple. 

3. Repeat step [2] until the Euclidean length of the segment (a„,6„) is smaller than e. 

4. Iterate the Poincare map for initial conditions a„ and 6„, obtaining a new pair (a„, 
Repeat this step until the length of the segment (a„, 6„) is larger than e. 

5. Repeat step [2j 

This method allows one to obtain a numerical trajectory which will stay very close to the 
chaotic saddle for an arbitrarily long period of time within an error e, and it is useful to 
compute statistical quantities which require long trajectories, such as the Lyapunov expo- 
nents. Its main disadvantage is that the PIM- Triple method may not work in the presence 
of two or more positive Lyapunov exponent^ (i.e., the stable manifold of the chaotic saddle 
has codimension higher than one). Instead, one can use other alternatives such as the PIM- 
simplex algorithmic, the Stagger-and-Step methocP^, or the gradient search algorithmP^l. 
Our implementation of the PIM-triple method follows the original description by Nusse and 
Yorke'^, and we set e = 10 



2. The sprinkler method 

In the sprinkler methocP^'^ the restraining region in the Poincare surface of section is 
covered by a two-dimensional (2D) grid of initial conditions, and then each initial condition 
is iterated until some time tc, which should be larger than the average escape time from the 
restraining region, and must be adjusted after some trial and error. The final points in the 
Poincare surface that remain in the restraining region after tc time units approximate the 
unstable manifold, their initial conditions approximate the stable manifold, and the points 



obtained at t = ^tc approximate the chaotic saddle. For most systems (including Eqs. (12) 



(|T5|)) ^ = 0.5 works fin(PS21. The accuracy of this method can be improved by choosing 



a large however, the number of initial conditions in the grid must be increased as well. 
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Our implementation of this method is based on the description given by Hsu et al.^, using 
a grid of 1000 x 1000 initial conditions, and we set tc = 20 time units. 

The chaotic saddle must coincide visually with the crossings of its stable and unstable 
manifolds obtained using the Sprinkler method. Since in our case the Poincare surface is 
a three-dimensional hypersurface, the visualization of the chaotic saddle and its stable and 
unstable manifolds can be affected by the choice of the 2D grid of initial conditions. We 
overcome this by choosing a grid of initial conditions in (oq, Oi), and setting the values of </> 
to = P{aQ), where P is a second-order polynomial which was obtained by plotting as a 
function of and computing a quadratic fit using the Poincare points of the chaotic saddle. 
This ensures that the chaotic saddle and its manifolds are entirely contained in the surface 
defined by (ao, ai, P(ao)), and the intersections between the stable and unstable manifolds 
coincide with the chaotic saddle. Details on this procedure will be provided in a future 
workpSI. 



IV. NONLINEAR ANALYSIS 
A. Parameter space 

We begin our numerical analysis by fixing 6 = —6 following Refs. [TT] and [15], and 



inspecting the stability of asymptotic solutions of Eqs. ( 12 )-( 15 ) as a function of the damping 
rates ui and z/2, i.e., in a two-dimensional section of the parameter space. Figure [l|^a) shows 
the value of the maximum non-zero Lyapunov exponent Amax in parameter space. This 
figure was obtained as follows 

1. Define a region of the parameter space by setting starting and ending values for ui 
and z/2. 



2. Set initial conditions and solve Eqs. (12)-(15) until t = T to discard transient values. 



3. Obtain Amax by computing the Lyapunov spectrum as explained in Section III A 



4. After t = Tmax, enough to ensure the convergence of Amax, increase the value of 1^2 
by a small amount Az/2 and use the last values of the integrated trajectory as initial 



conditions of Eqs. (12)-(15), to ensure that the same attractor is being "followed" 



across the parameter space. Repeat step |3] 
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5. After reaching the ending value of i>2, reset z/2 to its starting value, increase z/i by a 
small amount Az/i, and repeat from step [2] 

In our case we set {1^1,^2) e [6.2,7] x [6.2,7], T = 3000, Tmax = 20000, and Az/i = 
Ah'2 = 10^^. Positive values of Amax are represented using a colour scale, whereas negative 
values are shown using a gray scale. The colour scale varies from yellow (smaller values) 
to red (larger values), whereas the gray scale varies from white (smaller values) to black 
(larger values). Figure [l](a) clearly displays a central area of chaotic dynamics (Amax > 0) 
surrounded by ordered dynamics (Amax < 0). Embedded in the chaotic area are regular 
structures, enclosing periodic dynamics, identified as "shrimps" with a characteristic shape 
consisting of a main body and four elongated "antennae" . Fig. [l]^b) shows a detailed view 
of these self-similar structures, revealing a complex network of paths, or "window streets"'251 
composed by the entangled antennae and their intersections interconnecting shrimps. These 
paths allow one to move through the parameter space, while keeping the dynamics periodic. 
It is clear that the rich structures found in the parameter space are symmetric with respect 



to the line z/2 = z/i, which corresponds to the dynamics of the reduced system (16)-(18). 



B. Bifurcation diagram 

Next, we construct a bifurcation diagram by fixing z/i = 6.75 following previous workJ^^^, 
and choosing z/2 as our control parameter. Figure shows an overview of the bifurcation 
diagram using the amplitude of the first daughter wave ai as a function of z/2. This figure 



was constructed as follows. Set z/2 = 6.2, and set the initial conditions. Solve Eqs. (12)-(15), 
discarding the first 1500 transient Poincare points, and save the next 300 points, plotting 
ai against z/2. Then, increase z/2 by a small step, use the last Poincare point as the initial 
condition, and solve again the evolution equations until z/2 = 7. The resulting bifurcation 
diagram shown in Fig. [2]^a) agrees with the chaotic and periodic regimes seen in Fig. [l]^a). 
From Fig. [2]^a) we observe the presence of a limit cycle at z/2 = 6.2, represented by a 
period- 1 attractor. This means that the nonlinear waves described by the model equations 
display periodic dynamics. By increasing the control parameter the period of the limit cycle 
is doubled at z/2 ~ 6.34. At z/2 ~ 6.45 there is a sudden transition to a chaotic regime, which 
is occasionally interrupted by periodic windows as z/2 is increased. The largest periodic 
window occurs at z/2 ~ 6.5, and is indicated by an arrow. The bifurcation diagram becomes 
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very narrow with increasing z/2. At z/2 ~ 6.85 the solutions of the model equations converge 
asymptotically to a period-2 attractor. At U2 ~ 6.95 the periodicity of the attractor is 
reduced to 1. 

We compute the maximum non-zero Lyapunov exponent Amax for the same range of z/2. 
Fig. [2|^b) shows the value of Amax as a function of z/2. This figure is in agreement with 
Fig. [2|a) in the sense that periodic solutions correspond to negative values of Amax and the 
chaotic region is characterized by positive values of Amax- Periodic windows are reflected 
by sudden drops of the value of Amax towards negative values. Negative values of Amax for 
z/2 > 6.85 confirm that the system dynamics is periodic. 

We focus our analysis on the largest periodic window indicated by an arrow in Fig. |2| 
shown in detail in Fig. [3|^a). This periodic window is created by a tangent bifurcation at 
^2 = ^'SNB ~ 6.4845 known as a saddle-node bifurcation (SNB), where the chaotic attractor 
disappears and a stable and an unstable period-3 orbits are created. The stable periodic 
orbit (SPO) is displayed as a continuous line, while the unstable periodic orbit (UPO) is 
represented by a dashed line. The SPO undergoes a period-doubling cascade with increasing 
z/2, becoming a period-6 SPO at z/2 ~ 6.495, and a period-12 SPO at z/2 ~ 6.5025. Eventually, 
the periodicity of the SPO goes to infinity, becoming a chaotic attractor which can be 
distinguished as a set of three bands in Fig. [sj^a). At z/2 = z/jc ~ 6.5074 the chaotic 
attractor collides with the UPO created at the SNB and suddenly increases in size. This 
global bifurcation is known as an interior crises (IC), and marks the end of the periodic 
window. The UPO which mediates the IC is called a mediating UPO. 

The computation of the maximum Lyapunov exponent shown in Fig. [3|^b) clearly demon- 
strates the transition from chaotic to ordered states at SNB. At IC, there is a sudden increase 
in the value of Amax, for this reason the pre-crisis chaotic attractor with smaller Amax is 
called a weak chaotic attractor (WCA), while the enlarged, post-crisis chaotic attractor with 
larger Amax is called a strong chaotic attractor (SCA). 

After colliding with the mediating UPO the WCA loses stability, becoming a nonattract- 
ing chaotic set confined to the same banded region. For this reason it is called a banded 
chaotic saddle (BCS). Fig. [3|c) shows the same bifurcation diagram as in Fig. |3](a), from 
z/2 = 6.48 to z/2 = Vie- After IC, we remove the SCA and plot the evolution of the BCS as a 
function of z/2. The BCS was obtained using the PIM-triple procedure, and the restraining 
region is defined as the same region previously occupied by the WCA, bounded by the me- 
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diating UPO and its stable manifold. Fig. [sj^c) clearly shows the presence of empty regions 
or "gaps" that widen as 1^2 increases. 

Another type of chaotic saddle, which exists for the entire bifurcation diagram, is shown 
in Fig. ^d). This chaotic saddle is called a surrounding chaotic saddle (SCS), because it 
"surrounds" the space occupied by the attractors within the periodic window, and the space 
occupied by the BCS after the interior crisis. The SCS was obtained using the PIM-triple 
procedure. The restraining region was defined in different ways, depending on the value 
of z/2. For z/sNB < 1^2 < t^ic (i-e., inside the periodic window), in the periodic regime, we 
cover the SPO with a small disc of radius r = 10~^, and define the restraining region as the 
region of the phase space surrounding the SPO minus the disk. In the chaotic regime, the 
restraining region is defined as the region outside a rectangle that covers one of the bands of 
the WCA. For 1/2 > vic (i-e., after the interior crisis) we define the restraining region as the 
region located outside that occupied by the BCS. For z/2 < vs^b (i-e., before SNB) we use 
the same region as the first value of V2 in the periodic regime. Similarly to the BCS shown 
in Fig. [Sj^c), the SCS displays "gaps" that become wider with increasing 1/2, as displayed in 
Fig. gd). 

The SNB, which marks the beginning of the periodic window shown in Fig. |3| is charac- 
terized by a sudden transition from chaotic to periodic dynamics. Before SNB, the chaotic 
time series displays type-I Pomeau-Manneville intermittencjJ^ in which intervals resembling 
a periodic behavior are interrupted by occasional chaotic bursts. Figure |4|^a) shows the time 
series of the first daughter wave ai, represented by Poincare points, for U2 = 6.4844 < usnb, 
which clearly displays intermittent behavior. After SNB, the asymptotic dynamics con- 
verges to a period-3 solution, as evidenced by the time series shown in Fig. |4](b) for 
1^2 = 6.4845 > i^sNB- Recall that the period-3 stable periodic orbit created at SNB un- 
dergoes a period-doubling cascade with increasing z/2, and that the periodicity of the SPO 
eventually tends to infinity, leading to the WCA. In this regime the time series of Poincare 
points of the first daughter wave ai is characterized by weakly chaotic dynamics, and the 
points remain confined to three narrow bands as shown in Fig. |4]^c) for z/2 = 6.5074 < z/jc. 
After colliding with the mediating UPO at IC, which marks the end of the periodic window, 
the attractor increases in size and becomes the SCA. The time series of the daughter wave 
Oi now displays strong "bursts" which occasionally interrupt "laminar" dynamics that keep 
the memory of the weak chaotic regime, as observed in Fig. |4|^d) for 1^2 = 6.5075 > z/jc- 
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This intermittency is known as crisis-induced intermittenc}^^^'^. 

The chaotic bursts observed in the time series of Figure |4]^a), due to the type-I Pomeau- 
Manneville intermittency, become more frequent as the value of the control parameter de- 
parts from the SNB point, within the chaotic regime. This means that the average time r of 
the intervals resembling periodic behavior decreases with decreasing z/2. Figure |5](a) shows 
r as a function of the departure of z/2 from z/snb (black squares) in log-log scale. This figure 
confirms that r decreases as the "distance" between 1/2 and z/snb increases. Clearly, the av- 
erage time follows a power-law, with a slope 7 = —0.36 obtained by a linear fit represented 
by a continuous line in Fig. [5]^a). Similarly, the average time r in which the trajectory 
is confined to three bands as observed in the intermittent time series of Fig. ^d) due to 
the crisis-induced intermittency decreases with the departure of z/2 from z/jc in the chaotic 
regimd^, and follows a power-law behavior with a slope 7 = —0.62, as shown in Fig. IHFb). 



C. Saddle-node bifurcation 

The SNB marks the beginning of the period-3 window shown in Fig. [sj^a), where the 
chaotic attractor disappears and a stable periodic orbit and an unstable periodic orbit are 
created. In the Poincare plane, periodic orbits are represented by fixed points. The creation 
of fixed points across SNB can be elucidated by constructing return maps. Figure [6]^a) shows 
the value of the amplitude of the first daughter wave ai at the n-th intersection with the 
Poincare map, denoted by (ai)„, plotted against its third crossing with the map, represented 
by (ai)n+35 using points from a chaotic trajectory at z/2 = 6.48 < z/snb (thick line). Fixed 
points of periods 1 and 3 can be identified as intersections between the return map and 
the line given by (ai)„+3 = (ai)„ (dashed line). From Fig. [6]^a) the return map intersects 
the line at (ai)„ ~ 0.9923. The return map also approaches the line at (ai)„ ~ 0.9912, 
(oi)n ~ 0.9920, and at (ai)n ~ 0.99265. By increasing the value of the control parameter 
to z/2 = 6.4844 < z/sNB the return map becomes tangent to the line at the same three points 
(Fig. [6]^b)). For z/2 = 6.49 > z/snb, since the dynamics of the attractor changed from chaotic 
to periodic we construct the return map using points from the surrounding chaotic saddle 
(SCS, black points). Due to its fractal nature the SCS has empty spaces or "gaps" which 
prevent us from demonstrating the crossings between the return map and the dashed line. 
We solve this by fitting the Poincare points of the SCS to a curve using a cubic fit, and 
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superposing the interpolated points against the return map. The return map obtained by 
this technique is shown as a thin hne in Fig. ^c). It is evident that the interpolated points 
fit well the return map of the SCS, and the intersections between the return map and the 
line (ai)„ = (ai)„+3 clearly demonstrate the creation of a pair of stable-unstable period-3 
fixed points. 



D. Interior crisis 



The interior crisis is characterized by the collision between the WCA and the mediating 
UPO and its stable manifold at ICP^. A common technique to compute the stable manifold 
of an UPO requires the integration of the model equations backwards in timd^. However, 



Equations (12)-(15) turned out to be numerically unstable to be solved backwards in time 
with the traditional numerical integration schemes. Since the stable manifold of the mediat- 
ing UPO and the stable manifold of the chaotic saddle are practically indistinguishable^'^, 
an alternative approach is to compute the stable manifold of the SCS and extract its "con- 
tour" to represent the stable manifold of the mediating UPO near the collision site. Figure 
7| shows an overview of the stable (gray) and unstable (red) manifolds of the SCS (green), 
for z/2 = 6.5074 < uic- The WCA is shown in black, and the period-3 mediating UPO is 
represented by crosses. The SCS and its stable and unstable manifolds are obtained using 
the sprinkler method, in which a grid of 1000 x 1000 initial conditions in the (oq, Oi) plane 



was defined, and t = 10 as indicated in Section IIIB We use the following second-degree 



polynomial obtained by a quadratic fit to set the value of (p for the initial conditions 

= p(ao) = 5.315 - 0.44426ao + 0.029593a^. 



The chaotic saddles depicted in Figs. |3](c) and (d) have "gaps" clearly distinguished as empty 
spaces that widen with increasing z/2. After IC, these gaps are densely filled by the sudden 
appearance of a special set of UPOs termed coupling UPOs, because they have components 
in the surrounding and the band regions of the phase space, coupling themP^'SlESl^ This gap 
filling is an example of a bifurcation phenomenon called "explosion" whereby new recurrent 
points suddenly appear at a nonzero distance from any pre-existing recurrent point^. The 
gap filling is directly responsible for the crisis- induced intermittency observed in Fig. |4](d), 
because trajectories can "escape" from the banded region to the surrounding region via the 
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unstable manifold of the newly created coupling orbits. In Fig. |8[a) we show the post-crisis 
strong chaotic attractor (SCA, black) superposed by the period-3 mediating UPO (crosses) 
and its stable manifold (dashed line) obtained using the contour of the stable manifold of the 
SCS, for V2 = 6.515 > vic. In Fig. [8|b) we plot the BCS (black circles) and the SCS (gray 
circles) which constitute the so called geometrical "skeleton" of the post-crisis attractoi'^. 
Fig. [sj^c) shows one example of a period-10 coupling UPO (crosses). From this figure, it 
is clear that the coupling UPO has Poincare points located in both the surrounding region 
and the band region. This coupling UPO was detected by applying the Newton method to a 
large set of candidate points, obtained from the Poincare points of the SCA located within 
the gaps of either BCS or SCS. 

V. CONCLUSION 

In this paper we performed numerical simulations of a model of nonlinear modulational 
processes in a plasma, which describes the onset of Langmuir turbulence. We demonstrated 
the existence of shrimp-shaped structures in a two-parameter space. To our knowledge, this 
is the first time that these structures are detected in a model of nonlinear wave-wave inter- 
actions, rendering support to the universal nature of this phenomenon. Using two numerical 
methods, the PIM-triple method and the sprinkler method, we identified chaotic saddles and 
showed their evolution in the bifurcation diagram. We focused on two types of intermittency, 
the type-I Pomeau-Manneville intermittency, related to a saddle-node bifurcation, and the 
crisis-induced intermittency, related to an interior crisis. Both types of intermittency were 
shown to display power-law behavior as the control parameter departs from the transition 
point. After crisis, we were able to detect one example of a coupling unstable periodic orbit 
connecting two chaotic saddles, directly responsible for the intermittency observed in the 
time series after crisis. 

Nonlinear modulational processes can describe the formation of modulated wave en- 
velopes in plasmas and fluids which are often observed as intermittent "bursts" in a time 
series. For example, in astrophysical plasmas, observational data of radio signals from 
pulsars display strong bursts that are intermittent in timd^SESl Theoretical models of mod- 
ulational processes can also be used to explain the formation of extreme or "freak" waves in 
oceans, which are isolated, anomalously high-amplitude ocean waves, exceeding at least 2.2 
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times the significant wave height, which is defined as the average of the highest one-third of 
nearby waved^*^, and can appear intermittently in space or time, generally under rough sea 
condition d^^ l ^^ l In the open ocean, freak waves can be produced by nonlinear self-modulation 
of a slowly-varying wave trairP. Simple nonlinear solutions of the NLS equation can serve 
as prototypes of freak waves in the open ocean, as shown by several theoretical work^^^*^ 
and recent laboratory experiment^. A system of two coupled NLS equations were derived 
by Onorato et al.^, to describe the modulational instability in the ocean. Their model 
can explain the occurrence of freak waves in crossing sea states. Voronovich et al.^^ de- 
rived a NLS-type equation in which the bottom stress is taken into account, and found that 
the bottom friction can affect the formation of freak waves via a modulational instability. 
However, the three travelling waves truncation of Eq. ([T]) imposes some constraints on the 
applicability of this model to interpret the intermittent features observed in plasmas and 
fiuids. The choice of one pump wave and one pair of daughter waves might be justified in a 
weakly nonlinear regime, in which other sets of waves are strongly damped. In this regime, 
low- dimensional models of nonlinear wave-wave interactions can be used as an alternative 
to describe modulational processes in plasmas and fiuids. For example, Sanchez- Arriaga et 
al.'^ compared numerical simulations of the derivative nonlinear Schrodinger equation with 
simulations of a reduced model consisting of three travelling waves in a plasma, and found 
that the reduced model is in agreement with the DNLS when the energy in the system is 
small and right-hand polarized waves are present. Janssenl^ showed that, when the wave 
steepening is sufficiently weak, the theory of nonlinear four-wave interactions is in good 
agreement with direct simulations of an ensemble of ocean waves using Zakharov equations 
and nonlinear Schrodinger equation, and can explain the formation of freak waves. Hence, 
the model used in this paper may be useful to describe qualitatively the dynamics of non- 
linear modulational interactions. In summary, our work offers a theoretical explanation for 
the generation of intermittent bursts resulting from the interaction of chaotic saddles and 
coupling unstable periodic orbits, and provides new insights to the understanding of these 
nonlinear phenomena commonly observed in plasmas and fluids. 



18 



ACKNOWLEDGMENTS 



This work is supported by CNPq and FAPESP. A.C.L.C. is grateful for the award of 
a Marie Curie International Incoming Fellowship and the hospitality of Paris Observatory. 
The authors thank the referee for valuable comments. 



REFERENCES 

^H. Kim, R. L. Stcnzcl and A. Y. Wong, "Development of 'Cavitons' and trapping of rf 
field", Phys. Rev. Lett. 33, 15, 886-889 (1974). 

^F. Amiranoff, M. Laberge, J. R. Marques, F. Mouhn, E. Fabre, B. Cros, G. Matthieussent, 
P. Benkheiri, F. Jacquet, J. Meyer, Ph. Mine, C. Stenz, and P. Mora, "Observation of 
modulational instability in Nd-laser beat-wave experiments" , Phys, Rev. Lett. 68, 3710- 

3713 (1992). 

^R. E. Ergun, C. W. Carlson, J. P. McFadden, J. H. Clemmons, "Evidence of a transverse 
Langmuir modulational instability in a space plasma", Geophys. Res. Lett. 18, 7, 1177- 
1180 (1991). 

^J. Bonnell, P. Kintner, J.-E. Wahlund and J. A. Holtet, "Modulated Langmuir waves: 
observations from Freja and SCIFER", J. Geophys. Res. 102, A8, 17233-17240 (1997). 

^J. C. Weatherall, "Modulational instability, mode conversion, and radio emission in the 
magnetized pair plasma of pulsars", Astrophys. J. 483, 402-413, (1997). 

^P. A. E. M. Janssen, "Nonlinear four-wave interactions and freak waves", J. Phys. Ocean. 
33, 863-884 (2003). 

^M. V. Goldman, "Strong turbulence of plasma waves". Rev. Mod. Phys. 56, 709-735 

(1984). 

*A. Pataraya and G. Melikidze, "Modulational instability of nonlinear waves in the rcla- 
tivistic plasma with account of the nonlinear Landau damping" Astrophys. Space Sci. 68, 
61-71 (1980). 

^P. A. Robinson, I. H. Cairns, and N. I. Smith, "Unified theory of monochromatic and 
broadband modulational and decay instabilities of Langmuir waves", Phys. Plasmas 9, 10, 
4149 (2002). 

19 



^°F. A. Asenjo, F. A. Borotto, A. C.-L. Chian, V. Munoz, J. A. Valdivia and E. L. Rempel, 

"Self-modulation of nonlinear waves in a weakly magnetized relativistic electron-positron 

plasma with temperature", Phys. Rev. E 85, 046406 (2012). 
^^D. A. Russell and E. Ott, "Chaotic (strange) and periodic behavior in instability saturation 

by the oscillating two-stream instability", Phys. Fluids 24, 11, 1976-1988 (1981). 
^^S. Gosh, and K. Papadopoulos, "The onset of Alfvenic turbulence", Phys. Fluids 30, 

1371-1387 (1987). 

^^D. F. DuBois, H. A. Rose and D. R. Nicholson, "Statistical theories of Langmuir turbulence 

II: subsonic to sonic transition", Phys. Fluids 28, 202-218 (1985). 
"^^A. C.-L. Chian and J. R. Abaldc, "Fundamental plasma radiation generated by a travelling 

Langmuir wave: hybrid stimulated modulational instability" , J. Plasma Phys. 57, 4, 753- 

763 (1997). 

^^R. A. Miranda, E. L. Rempel, A. C.-L. Chian, and F. A. Borotto, "Intermittent chaos 
in nonlinear wave-wave interactions in space plasmas", J. Atmos. Sol.-Terr. Phys. 67, 
1852-1858 (2005). 

^^J. R. Sanmartin, O. Lopez-ReboUal, E. del Rio and S. Elaskar, "Hard transition to chaotic 

dynamics in Alfven wave fronts", Phys. Plasmas 11, 2026 (2004). 
^'^G. Sanchez-Arriaga, J. R. Sanmartin, and S. Elaskar, "Damping models in the truncated 

derivative nonlinear Schrodinger equation", Phys. Plasmas 14, 082108 (2007). 
^^G. Sanchez-Arriaga, T. Hada, and Y. Nariyuki, "The truncation model of the derivative 

nonlinear Schrodinger equation", Phys. Plasmas 16, 042302 (2009a). 
^^G. Sanchez-Arriaga, T. Hada, and Y. Nariyuki, "Truncation model in the triple-degenerate 

derivative nonlinear Schrodinger equation" Phys. Plasmas 16, 042303 (2009b). 
^"N. Mori and P. A. E. M. Janssen, "On kurtosis and occurrence probability of freak waves" 

J. Phys. Ocean. 36, 1471-1483 (2006). 
^^N. Mori, M. Onorato, P. A. E. M. Janssen, A. R. Osborne and M. Serio, "On the extreme 

statistics of long-crested deep water waves: Theory and experiments", J. Geophys. Res. 

112, C09011 (2007). 

^^J. A. C. Gallas, "Structure of the parameter space of the Henon map", Phys. Rev. Lett. 
70, 2714-2717 (1993). 

^^Y. Zou, R. V. Donner, J. F. Donges, N. Marwan and J. Kurths, "Identifying complex 
periodic windows in continuous-time dynamical systems using recurrence-based methods" , 



20 



CHAOS 20, 043130 (2010). 

R. Stoop, P. Benner, and Y. Uwate, "Real-world existence and origins of the spiral orga- 
nization of shrimp-shaped domains", Phys. Rev. Lett. 105, 074102 (2010). 
'E. N. Lorenz, "Compound windows of the Henon map", Physica D 237, 1689-1704 (2008). 
'C. Bonatto, J. C. Garrcau, J. A. C. Gallas, "Self-similarities in the frequency-amplitude 
space of a loss-modulated CO2 laser", Phys. Rev. Lett. 95, 143905 (2005). 

D. M. Maranhao, M. S. Baptista, J. C. Sartorelli, and I. L. Caldas, "Experimental obser- 
vation of a complex periodic window", Phys. Rev. E 77, 037202 (2008). 

'C. Grebogi, E. Ott, J. A. Yorke "Crises, sudden changes in chaotic attractors, and transient 
chaos", Physica D, 7 , 181-200 (1983). 

'H. Kantz, P. Grassberger, "Repellers, semi-attractors, and long-lived chaotic transients" , 
Physica D 17 75-86 (1985). 

'G.-H. Hsu, E. Ott, C. Grebogi, "Strange saddles and the dimensions of their invariant 
manifolds", Phys. Lett. A 127 4, 199-204, (1988). 

Y.-C. Lai and T. Tel, Transient Chaos: Complex Dynamics on Finite-Time Scales, 
(Springer, 2011). 

'A. C.-L. Chian, R. A. Miranda, E. L. Rempel, Y. Saiki, M. Yamada, "Amplitude-phase 
synchronization at the onset of permanent spatiotemporal chaos" , Phys. Rev. Lett. 104 , 
254102 (2010). 

'E. L. Rempel, A. C-L. Chian, E. E. N. Macau, R. R. Rosa, "Analysis of chaotic saddles 
in low-dimensional dynamical systems: the derivative nonlinear Schrodinger equation", 
Physica D 199 , 407-424 (2004). 

E. L. Rempel, A. C.-L. Chian, "Origin of transient and intermittent dynamics in spa- 
tiotemporal chaotic systems", Phys. Rev. Lett. 98 , 014101 (2007). 

'D. R. Nicholson, Introduction to plasma theory, (John Wiley and Sons, New York, 1983). 
'F. F. Chen, Introduction to plasma physics and controlled fusion, (Plenum Press, New 
York, 1984). 

L. R. Petzold, "Automatic selection of methods for solving stiff and nonstiff systems of 
ordinary differential equations", Siam J. Sci. Stat. Comput. 4, 136 (1983). 
'V. I. Oseledec, "A multiplicative ergodic theorem, Lyapunov characteristic numbers for 
dynamical systems". Trans. Moscow Math. Soc. 19, 197-231 (1968). 



21 



^^I. Shimada and T. Nagashima, "A numerical approach to ergodic problem of dissipative 
dynamical systems", Prog. Theor. Phys. 61, 1605-1616 (1979). 

^°G. Benettin, L. Galgani, A. Giorgilli and J.-M. Strelcyn, "Lyapunov characteristic expo- 
nent for smooth dynamical systems and for Hamiltonian systems; A method for computing 
aU of them", Meccanica 15, 9-20 (1980). 

^^W.-H. Steeb, Linearization procedure and nonlinear systems of differential and difference 
equations, in Nonlinear Phenomena in Chemical Dynamics, edited by C. Vidal and A. 
Pacault (Springer, Berhn, 1980), p. 275. 

^^W.-H. Steeb and J. A. Louw, Chaos and quantum chaos (World Scientific, Philadelphia, 
1986). 

^^E. Ott, Chaos in dynamical systems (Cambridge Univ. Press, New York, 1993). 

^'^A. Wolf, J. B. Swift, H. L. Swinney and J. A. Vastano, "Determining Lyapunov exponents 
from a time series", Physica D 16, 285-317 (1985). 

^^H. E. Nusse, J. A. Yorke, "A procedure for finding numerical trajectories on chaotic sad- 
dles", Physica D 36 137-156 (1989). 

^^P. Moresco, S. P. Dawson, "The PIM-simplex method: an extension of the PIM-triple 
method to saddles with an arbitrary number of expanding directions". Physica D 126 
38-48 (1999). 

^'^D. Sweet, H. E. Nusse, J. A. Yorke, "Stagger-and-Step method: detecting and computing 

chaotic saddles in higher dimensions". Phys. Rev. Lett. 86 2261-2264 (2001). 
^^E. M. BoUt, "The path towards a longer life: on invariant sets and the escape time 

landscape". Int. J. Bif. Chaos 15 1615-1624 (2005). 
^^R. A. Miranda, E. L. Rempel and A. C.-L. Chian, "Chaotic saddles in nonlinear three-wave 

interactions in a plasma" , Physica D, to be submitted. 
^'^Y. Pomeau, P. Manneville, "Intermittent transition to turbulence in dissipative dynamical 

systems". Comm. Math. Phys. 74 189-197 (1980). 
^^C. Grebogi, E. Ott, F. Romeiras, J. A. Yorke "Critical exponents for crisis-induced inter- 

mittency", Phys. Rev. A, 36 , 11, 5365-5380 (1987). 
^^Z. You, E. J. Kostelich, J. A. Yorke, "Calculating stable and unstable manifolds". Int. J. 

Bif. Chaos 1, 605-623 (1991). 
^^E. M. Ziemniak, C. Jung, T. Tel, "Tracer dynamics in open hydrodynamical flows as 

chaotic scattering", Physica D 76, 123-146 (1994). 



22 



^ K. G. Szabo, Y.-C. Lai, T. Tel, C. Grebogi, "Topological scaling and gap filling at crisis", 

Phys. Rev. E 61, 5019-5032 (2000). 
^^K. G. Szabo, Y.-C. Lai, T. Tel, C. Grebogi, "Critical exponent for gap filling at crisis", 

Phys. Rev. Lett. 77, 3102-3105 (1996). 
^^C. Robert, K. T. AUigood, E. Ott, J. A. Yorke, "Explosions of chaotic sets", Physica D 

144, 44-61 (2000). 

^''K. G. Szabo, T. Tel, "Transient Chaos as the backbone of dynamics on strange attractors" , 

Phys. Lett. A 196 , 173-180 (1994). 
^^J. S. Deneva et al., "Arecibo pulsar survey using ALFA: probing radio pulsar intermittency 

and transients", Astrophys. J. 703, 2259-2274 (2009). 
^^A. C.-L. Chian, C. F. Kennel, "Self-modulational formation of pulsar microstructures" , 

Astrophys. Space Sci. 97, 9-18 (1983). 
^°N. Mori, P. C. Liu and T. Yasuda, "Analysis of freak wave measurements in the Sea of 

Japan", Ocean Engineering 29 , 1399-1414, (2002). 
^^A. R. Osborne, M. Onorato, M. Serio, "The nonlinear dynamics of rogue waves and holes 

in deep-water gravity wave trains", Phys. Lett. A 275, 386-393 (2000). 
^^Andonowati, N. Karjanto, E. van Groesen, "Extreme wave phenomena in down-stream 

running modulated waves", App. Math. Model. 31, 1425-1443 (2007). 
^^D. H. Peregrine, "Water waves, nonhnear Schrodinger equation and their solutions", J. 

Aust. Math. Soc. Series B, Appl. Math. 25, 16-43 (1983). 
^^V. V. Voronovich, V. I. Shrira, and G. Thomas, "Can bottom friction suppress 'freak wave' 

formation?", J. Fluid Mech. 604, 263-296 (2008). 
^^A. Chabchoub, N. P. Hoffmann and N. Akhmediev, "Rogue wave observation in a water 

wave tank", Phys. Rev. Lett. 106, 204502 (2011). 
^^M. Onorato, A. R. Osborne and M. Serio, "Modulational instabihty in crossing sea states: 

a possible mechanism for the formation of freak waves", Phys. Rev. Lett. 96, 014503 

(2006). 



23 




FIG. 1. (Color online) Shrimp-shaped self-similar structures in a two-parameter space. Upper 
panel: Maximum non-zero Lyapunov exponent Amax as a function of ui and 1^2- The colour scale 
indicates chaotic dynamics, whereas the gray scale denotes periodic oscillations. Lower panel: 
Enlargement of a subregion of (a). 
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Figure 2 
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FIG. 2. (a) Bifurcation diagram for ai as a function of f2. The arrow indicates a period-3 window, 
(b) Maximum non-zero Lyapunov exponent Amax as a function of V2. 
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FIG. 3. (a) Bifurcation diagram for ai as a function of 1/2. The dashed hnes represent the evolution 
of the period-3 unstable periodic orbit, SNB denotes a saddle-node bifurcation and IC denotes 
an interior crisis, (b) Maximum non-zero Lyapunov exponent Amax as a function of ^2- (c) 

Bifurcation diagram showing the conversion of the weak chaotic attractor (black) to a banded 
chaotic saddle (gray) and its evolution after crisis, (d) Bifurcation diagram showing the evolution 
of the surrounding chaotic saddle. 
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Figure 4 
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FIG. 4. (a) Time series of the Poincare points of ai at V2 = 6.4844, showing the characteristics of 
a type-I Pomeau-Manneville intermittency. (b) Time series of the Poincare points of ai at V2 = 
6.4845 > I'sNB, showing period-3 dynamics, (c) Time series of pre-crisis weak chaos represented 
by the Poincare points of ai at 1^2 = 6.5074 < i/jc- (d) Time series of ai at 1^2 = 6.5075 > fjc 
represented by the Poincare points, showing a crisis-induced intermittency alternating episodically 
between periods of weak and strong chaos. 
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FIG. 5. (a) Average time r between intermittent events (black squares) as a function of the 
departure of 1^2 from i^snb in the chaotic regime before the saddle-node bifurcation. The straight 
line represents a linear fit with a slope 7 = —0.36. (b) Average time r between intermittent bursts 
(black squares) as a function of the departure of 1^2 from V[c in the chaotic regime after the interior 
crisis, and a linear fit (straight line) with a slope 7 = —0.62, in a log-log scale. 
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Figure 6 
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FIG. 6. Return map ((ai)„, (ai)„+3) for (a) a chaotic attractor at 1/2 = 6.48 < usnb- The dashed 
hne represents (ai)„ = (ai)„+3. (b) Chaotic attractor at U2 = 6.4844 < vsnb- (c) Surrounding 
chaotic saddle at U2 = 6.49 > usnb (black dots) and a curve obtained using a cubic fit (thin line). 
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Figure 7 




FIG. 7. (Color online) Poincare plots of the surrounding chaotic saddle (green) right before crisis 
(z^2 = 6.5074 < i/ic) and its stable (gray) and unstable (red) manifolds. The weak chaotic attractor 
is represented by black thick lines, and the period-3 mediating unstable periodic orbit by blue 
crosses. 
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Figure 8 
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FIG. 8. Coexistence of a chaotic attractor, chaotic saddles and a coupling orbit after the interior 
crisis (f2 = 6.515 > i^c)- (a) Strong chaotic attractor (SCA) and the mediating unstable periodic 
orbit (UPO) represented by crosses. The dashed lines represent the stable manifold of the mediating 
UPO. (b) Banded chaotic saddle (BCS, black) and surrounding chaotic saddle (SCS, gray), (c) 
BCS and SCS superposed by a period-14 coupling unstable periodic orbit. The banded region is 
denoted by B and the surrounding region by S. 
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